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Abstract 

In this paper we compare dealing methods and we generahze the results of Conger and Howald. They 
compared three well known dealing methods like the ordered dealing, the cyclic dealing and the back and 
forth dealing after an a shuffle in the case when there are 52 distinct cards in the deck and four players. 
In this paper we analized how fast the different dealing methods randomize the deck in the case when 
there are repeating cards. The main generalization of our work is that we compared these methods for 
those cases when the deck contains repeated cards of an arbitrary number of different types and only 
the hands dealt to players are of interest, moreover the number of cards in the deck and the number of 
players are also arbitrary. For the problem we used the variation distance that we stripped in Taylor 
series in the variable . We do that for different dealing methods, and for an arbitrary number of types 
of cards, number of players, and size of the deck. Denoting by s the number of cards dealt to each player, 
we prove for an arbitrary initial deck that the leading term (i.e., the coefficient of a~^) differs by a factor 
of s when comparing the ordered and the cyclic dealing. For odd s values it again differs by a factor of 
s when comparing the cyclic and the back and forth method, while this leading term is zero in the back 
and forth method for even s values. These differences mean that if a large number of riffie shuffies are 
performed, changing to the cyclic method instead of the ordered one (or to the back and forth method 
instead of the ordered one if s is odd) saves approximately logj s riffie shuffles to achieve the same level 
of randomness. We also show numerical values of the leading term in a couple of explicit scenarios. 



1 Introduction 

In card-games a very important requirement is that, after the mixing procedure, every hand dealt to players 
should have approximately the same probability. The main question is: "how fast can the riffling and shuffling 
method randomize the deck?" To examine this question, we need some definitions. 

In 1955 Gilbert and Shannon [S] used the riffle shuffling by examining the shuffling method. In the 1980s 
Reeds [B] and Aldous [T] added the assumption that every possible cut / rifEe combinations have equal chances 
using the riffle shuffling, and that has become known as the Gilbert-Shannon Reeds or GSR model of card 
shuffling. Riffle shuffling is the most fashionable way to randomize the initial deck. First we have to cut the 
deck into two packets. If the first packet contains k and the second packet n — k cards then the probability 

of this cut is exactly . After the cut we have to mix the packets together such that the cards of each 
packet maintain their relative order. We suppose that each possible interleaving is equally likely. There are 
(2) possible interleavings (we only have to choose the k places among n places for the first packet, then 
the other n — k places are determined). So we get that any possible shuffling where k is the size of the cut 

has ^ • — probability. This value does not depend on fc, so each pair of a cut and interleaving is 

equally likely. In 1992, Bayer and Diaconis |2] analized a generalization of the riffle shuffling. They introduce 
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the a-shuffle to make the mixing problem easier. First we cut the deck into a packets. If the sizes of the 

packets are Pi,P2, ■ ■ ■ ,Pk then the probabiHty of this particular cut is ipiiPiipiP^?! . After this cut we mix the 
cards together such that the cards of each packet maintain their relative order. If we imagine an a-handed 
dealer, who can interleave the cards like a two-handed dealer, we can see what a-shuffle means. They prove 
that making a randomly selected a-shuffle and then making a randomly selected &-shuffle is equivalent to 
making a randomly selected a ■ b shuffle, if we add the assumption as before that every possible cut/riffle 
combinations are equally likely. In particular, this implies that a sequence of i riffle shuffles is equivalent to 
a single 2*-shuffle. Bayer and Diaconis used the variation distance to analyze the card mixing problem. This 
distance is defined as: 

\K-U\\:^\ ^ K{n)-Uin)\, 

where Pa(7r) is the probability of a particular permutation tt after an a-shuffle, Sn is the symmetric group 
of degree n, U represents the uniform distribution on permutations (J7(7r)=^ for all tt S Sn), each card is 
distinct. Initially, the deck is ordered: if we label the cards by 1, 2, . . . , n, then in the initial deck card i 
will be in position i. Bayer and Diaconis found an explicit formula for P^: 

^ ^ a" V n 

where n is the size of the deck and 



des(7r) := #{i : 7r(i) > Tr{i + 1)}. 

In our case we consider a permutation from 1,2, ...,n, to itself, so if we apply tt to a sequence of objects, 
then the object in position i goes to position 7r(i). For example, the permutation 7ri=[43125] changes our 
initial ordering to 34215, as well as rearranging 25431 to 43521, and 53412 to 41352. We can easily check 
that des(7ri)=2. An interesting generalization is when we allow two cards to have the same value. This 
makes our problem more complicated because decks (ordered sequence of cards) and transformations cannot 
be identified with permutations. On the other hand there is a set of permutations for each pair which 
transform the first into the second. Another problem is that the initial order of deck affects how fast the 
distribution approaches uniform. If D' is some rearrangement of D, let S{D,D') be the set of permutations 
which transform D into D' . Thus the transition probability between D and D' is 

Pa{D ^ D') :- P-W- 

7reS(D,D') 

If we consider our explicit formula, we get the following: 

™ /„ „/N 1 , /a + rt — d— 1\ 

d ^ ^ 

where is the number of permutations in 5(-D, D') with d descents. 

Conger and Viswanath [4 proved that the calculation of transition probabilities is a #P-complete prob- 
lem. Most people believe that #P-complete problems do not admit efficient solutions, so a possible way to 
examine this question is to approximate this probability. Conger and Howald [3^ proved a theorem, which 
allows us to approximate transition probablitites, when a is large, given some simple information about the 
decks. To understand the theorem, we have to make some definitions. 

If a and 6 are card values, we say that D has an a — & digraph at i if -D(i) = a and D(i + 1) — 6. We say 
that D has an a — 6 pair at (i, if i < j, = a, and D{j) = b. Let 

W{D, a, b) :— i^{a — b digraphs in D} — #{6 — a digraphs in D}, 

Z{D, a, b) :~ ^{a — b pairs in D} — =ff{b — a pairs in D}. 
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For example, the following deck, that consists of red (i?) and black (B) cards, has \ R~ B digraph, 2 B — R 
digraphs, 12 R — B pairs, 13 B — R pairs: 



D := BRRRBBBBRR. 

So W{D, B,R) = 2-1 = 1, Z{D, B, i?) = 13 - 12 = 1. 

We can easily see that W and Z are antisymmetric in a and b: 

W{D, a, b) = -W{D, 6, a), Z{D, a, b) = -Z{D, b, a). 

Conger and Howald [3^ proved the following theorem: 

(1) Va{D -^D') = ^+ c,{D,D')a-^ + 0{a-^), 
where N is the number of reorderings of D and 

(2) c,(i^.i^o^^y^(^'"'^)^(^>'^\ 

2N ^ UaTlb 
a<b 

where Ua is the number of cards of value a, rib is the number of cards of value b. 

In this paper we will analyze dealing methods in the case when there are repeated cards in the deck and 
only the hands dealt to players are of interest. Conger and Howald [3] proved that switching from ordered 
dealing to cyclic dealing improves the randomness by a factor of 13, and switching from cyclic dealing to 
back and forth dealing also improves the randomness by a factor of 13, when there are 4 players and 52 
distinct cards. In Section 2 we describe the case when there can be repeated cards in the deck. In Section 
3 we analyze arbitrary dealing methods with repeated cards, in section 4 we analyze those cases when there 
are more players. We will write up a formula which shows us the effectiveness of a dealing method in the 
case when only the hands dealt to players are of interest, and we will prove from this formula, that switching 
from ordered dealing to cyclic dealing improves the randomness by a factor of s, and switching from cyclic 
dealing to back and forth dealing improves the randomness by a factor of s, when s is an odd number, thus 
switching from cyclic dealing to back and forth dealing improves the randomness in a huge amount, when s 
is an even number. Here s is the number of cards each player receives. In section 5 we analyze cases when 
the initial deck is arbitrary and we prove the former results also for these cases. Section 6 considers explicit 
computations when there are only two or three types of cards. 



2 The basics of our model 

First we consider the following case. Our deck contains 4s cards and there are k different types of cards: 
^1) -F2) • • • ) Pk- In the initial deck the first pi cards are Pi coloured, the next p2 cards are P2 coloured, • • • , 
and the last pk cards are Pk coloured (X^iLiPi — 4s). We a-shufhe the deck then we deal it out to the four 
players. The four players are North(N), East(E), South(S) and West(W). We use the variation distance as 
a level of randomness of the deck: 

l|P<^-n||:=^EliP"H-nMI' 

where ft is the set of possible configurations. 

fl = {pi,P2, ■ ■ ■ ,Pk ■ Pi.N,Pi,E,Pi.,SiPi.w > 0'-,P2,N,P2.E,P2,s,P2,w > 0; . . . ; Pk,N , Pk,E, Pk.,s i Pk,w > 0; 

Pl,N +Pl,E +Pl,S +Pl,W = Pi;P2,N +P2,E + P2,S + P2,W = P2] ■ ■ ■;Pk,N + Pk,E + Pk,S + Pk,W = Pk] 
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'^Pt.N = s;^p.,^E = s;^pi^s = s;^p,^w = s}, 



i=l 



where ft^jv, Pi^s, Pi,s^ Pi,w is the number of Pi coloured cards which North, East, South, West has respectively 
after the dealing method, Pi = iPi,N,Pi,E,Pi,s,Pi,w)- 

n(w) represents the stationary disribution, which supposes that each permutation of the 4s cards is 
equally likely. 

k 



Pt 



(4s)! ^J^ Pi^N^. -pi^E^. -Pi^s^. -Pi^w 

\-n||:=i^|p,H-nHH 



2 ^ 



Va{D^ D')-Uiuj) 

D':{pi,p2,--- ,Pk)=uj 



Here we suppose that deck D is the initial sequence of cards, and D' is some rearrangement of D. 
For computing Pa(-D — !• 13') we use ([T|): 

¥a{D -^D') = ^+ c,{D,D')a-^ + 0{a-^). 



In our case N — 
deck. 

We easily get the following: 



(4s) 



We assume an implicit order on the values according to the inital 



(3) 



D':{pi,P2,--- ■,Pk)=l^ « — 1 



Pi,N^- ■Pi,E^- ■Pi,S^- ■Pi,W^- 



Using these connections we get: (see also (|2 



J2 Va{D^D') 

D':{pi,p2.,--- ..Pk)=i^ 



E ^+«"' E c,{D,D')+ 0{a-') 



N 



D':(pi,P2,--- ,Ph)=^ 
k \ 



D':{pi,P2,--- ,Pk)=t^ 
-1 

(4s)! 



E T^v n^.' E 



2^ I „-i '^^ 



\a<b 



D':(pi,p2,--- ,Pk)='^ 

W{D,a,b)Z{D',a,h)^ 
nam I 



n(w) + 0(a"^) +a 



(4s) 



■Qp , j I ^ Y.D':(p-^, p-2,- ,p-k)=uj ^(^'^ Pj+l) 
\3 = l ) \i=\ 



Pi ■ Pz+1 



because Wifi, Pi,P,+i) = 1 (i = 1, • • • , fc - 1), and W{D, A, B) = if B =^ {A + 1) by virtue of the initial 
deck. 

Thus the variation distance will be the following: 



(4) i Y - n(.)| = a-^^_[l[p,^] Y 



wen 



E 



Pi ■ Pi+i 



+ 0(a-2). 
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3 Dealing methods 



Now we consider the dealing methods. The dealer o-shuffles the deck and then deals it out to the four 
players. We can describe the dealing method function as a function / from 1, 2, ... 4s to A'', E, S, W. 
For example: 

/(I) = /(2) = /(3) = • • • = /(,s - 1) = f{s) = N, 
f{s + 1) = /(s + 2) - /(s + 3) = • • • = /(2s - 1) = /(2s) = E, 
/(2s + 1) - /(2s + 2) = /(2s + 3) = . . . = /(3s - 1) = /(3s) = S, 
/(3s + 1) = /(3s + 2) = /(3s + 3) = • • • = /(4s - 1) = /(4s) = W, 

represents "ordered dealing", where the first s cards go to North, the next s cards go to East, the next s 
cards go to South and the last s cards go to West. We can represent the dealing method with a sequence of 

N, E, S and W. In this case the representing sequence is: 

NNN . . . NN EEE ...EE SSS ...SS WWW . . . WW. 
Another dealing method is the cyclic dealing: 



/(I) 


= /(5) 


= /(9) = -- 


• = /(4s- 


7) = 


= /(4s- 


3) = N, 


/(2) 


= /(6) 


= /(lO) = • 


• • = /(4s 


-6) 


= /(4s 


-2) = E, 


,/(3) 


= /(7) 


-,/(ll) = - 


• • = /(4s 


-5) 


= /(4s 


-1)^S, 


/(4) 


= /(8) 


= /(12) = • 


• • = /(4s 


-4) 


= /(4s) 


^W, 



so the top card goes to North, the second goes to East, the third goes to West, etc. The representing sequence 
is: 

NESWNESWNESW . . . NESWNESW. 
The back and forth dealing is the following for even s numbers: 

/(I) = /(8) = /(9) = /(16) = /(17) = /(24) = /(25) = /(32) = • • • = /(4s) = N, 
/(2) = /(7) = /(lO) = /(15) = /(18) = /(23) = /(26) = /(31) = • • • = /(4s -1)=E, 
/(3) = /(6) = /(II) = /(14) = /(19) = /(22) = /(27) = /(30) = • • • = /(4s - 2) = 5, 
/(4) = /(5) = /(12) = /(13) = /(20) = /(21) = /(28) = /(29) = • • • = /(4s -^) = W. 

The representing sequence is the following in this case: 

NESWWSENNESWWSEN . . . NESWWSEN. 

The back and forth dealing is the following for odd s numbers: 

/(I) = /(8) = /(9) = /(16) = /(IT) = /(24) = /(25) = /(32) = ■ ■ • = /(4s - 3) = iV, 
/(2) = /(7) = /(lO) = /(15) = /(18) = /(23) = /(26) = /(31) = • • • = /(4s -2) = E, 
m = m = /(II) = /(14) = /(19) = /(22) = /(27) = /(30) = • • • = /(4s -1) = S, 
/(4) = /(5) = /(12) = /(13) = /(20) = /(21) = /(28) = /(29) = • • • = /(4s) = W. 

The representing sequence is the following in this case: 

NESWWSENNESWWSEN . . . NESW. 

We will write up a formula which allows us to compute the coefficient of for any deal method, if we 
know the dealing method function. Here we suppose that in the initial deck the first type is black (B), the 
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second type is red (R). In this case b is the number of black cards in the deck, r is the number of red cards 
in the deck and p is the number of non-red and non-black cards in the deck. 

Let Ni, El, Si, Wi be the first position which goes to North, East, South, West respecively in the dealing 
method. Let N2, E2, S2, W2 be the second position which goes to North East, South, West respectively in 
the dealing method, etc. Let pjv, Pe, PSi Pw be the number of non-red and non-black cards which North, 
East, South, West has respectively after the dealing method. Let 6jv, bs, bw be the number of black 
cards which North, East, South, West has respectively after the dealing method. 



Theorem 3.1. 



E 



Z{D',B,R) 

D':(b,r,p3,P4,,--- ,Pk)=uJ 



^f-JlPt,N^- ■ Pi,E^- ■ Pt,S^- ■ Pi,wl^ ^ SV S p^l ^ 

bE (PN ■ Z{N, E)+ps- Z{S, E)+pw Z{W, E) 



+ ■ 
s 



* p=i 



^ bs / PN ■ Z{N, S)+pE- Z{E, S)+pw Z{W, S) 
s \ s 

p=l 

^ bw ^ PN ■ Z{N, W)+pE- Z{E, W)+ps- Z{S, W) ^ ^ ^^^^^ ^^ 



where 



Z{hj) = #{i — j pairs in the representing sequence of the dealing method} 
— #{j — i pairs in the representing sequence of the dealing method}. 

Proof. First we will prove that for a particular permutation D': 
(5) 

4s 

Z{D', B,R) = ^{4s + l- 2j+(the number of non-red and non-black cards in D' before the i*^ position) 

— (the number of non-red and non-black cards in D' after the i^^ position))- 
l{in the i^^ position there is a black card in permutation D'}. 

If we change in the position i the value from red to black in a deck D', then we can compute the 
change of the value Z{D',B,R). Put the case that in the first i — 1 cards there are A black, C red cards 
and G non-red and non-black cards; the card in the position i is red coloured; and in the last 4s — i 
cards there are D black, F red cards and H non-red and non-black cards. SoA + C + G = i— 1 
and D + F + H — As — i. In this case the card in the position i stands in {B — R pairs} with those 
A cards in the first i — 1 cards which are black coloured; and stands in {R — B pairs} with those D 
cards in the last 4,s — i cards which are black coloured. If we change in the position i the vahie from 
red to black we get that the change of the value Z{D',B,R) is: -A - C + D + F = -(2-l)+G'-|- 
4s — i — H = 4s + 1 — 2?' + (the number of non-red and non-black cards in D' before the position) — 
(the number of non-red and non-black cards in D' after the z*'^ position). We note that we only have to 
count the change between the card in the position i and the remaining cards. The A {B — R pairs} and the 
D {R — B pairs} vanish and there will be C new {R — B pairs} and F new {B — R pairs}, so the change 
of the value Z(D' , B, R) is indeed: —A — C + D + F. Let the initial deck have all its cards red colored, so 
in the inital deck the value Z{D' , B, R) is 0. First we change one by one in those positions the value from 
red to another card value in the initial deck in which the card value is non-black and non-red in the deck. 
After these changes the value Z{D' ,B,R) will be still in this revised deck. After this procedure, we have 
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to change one by one in those positions the value from red to black in this revised deck in which the card 
value is black in the deck D' , and then we can count one by one the change of the value Z{D' , B,R). If we 
sum up these changes we get the value Z{D' ,B,R), because in the initial deck and in the revised deck the 
value Z{D\ B, R) is 0. So we proved 

J2 Z{D',B,R)^ (45 + 1)6 

D':(b,f,p3.p4,--- ,Pk)=uj D':(b,f,p3,p4.--- .,Pk)=ul 

4s 

+ ^^(J7(z))l{in the i*'' position there is a black card in permutation Z?'}, 

D':(b,f,p3,pi,--- ,pk)=l^ ^=1 

where U{i) = (the number of non-red and non-black cards in D' before the i'^ position— the number of 
non-red and non-black cards in D' after the i*^ position— 2i). 

In order to compute X]_d' (f^(*))i introduce an auxiliary probability measure. According to this measure 
each permutation is equally likely. With the help of this measure we handle the sum as a conditional 
expectation. In this probability measure U{i) is a function of the permutation and, as such, becomes a 
random variable in our auxiliary probability measure, so we write this sum as a conditional expectation. 

(C/(i))l{in the i"^ position there is a black card in permutation D'} 

D':(b,f,p3,p4,--- .,Pk)=l^ 

\D' : (6, f, jr3, jr4, • • • = a;, in the i^^ position there is a black card| 

•E(C/(i)|in the i"^ position there is a black card in permutation D' , (&, f, P3,p4, • • • ,pk) — w) 
for i = 1, . . . , 4s. 

The justification of the last equality is that for each permutation Z?', where (6, f,p3,p4, • • • ,pk) — uj 
realized we sum up the value of U{i), if position i is black in permutation D' . In 

E(t/(i)|in position i there is a black card in permutation D\ {b,f,p3,p4, ■ ■ ■ ,pk) — uj), 

we also sum up these terms, but we divide each term by 

\D' : (6, r,p3,p4, • • • ,pk) — uj, in position i there is a black card|. 

Using these results we get the following: 

Z{D',B,R) = \D' : (5, f,p-3,p-4, • • • - c^|(4s + 1)6 

D':(h,f,p3,pi,--- ,Pk)=l^ 
4s 

\D' : (6, f,p3,p4, • • • ,pk) — UJ, in the i"^ position there is a black card| 

1=1 

•E([/(i)|in the i"^ position there is a black card in permutation D' , {b,f,p3,p4^, ■ ■ ■ ,pk) — uj) 

If we know the dealing method, we also know which player position i belongs to. Suppose that position 
i belongs to player North and in the i*^ position there is a black card. In this case: 

k 

\D' : (6,f,p-3,p-4, ■■■ ,p-k)^Lu\^bN- {s - 1)! • TT ^ , . 
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Suppose that position i belongs to player East and in the i*^ position there is a black card. In this case: 

k 

\D' : {b,f,p3,P4, ■ ■ ■ ,p-k) =u:\^hE- {s - 1)! • .s!^ [] 



J-- 

th 



Suppose that position i belongs to player South and in the i position there is a black card. In this case 

fe ^ 
\D' : (5,f,p-3,p-4, • • • ,Pk) = = 6s • (s - 1)! • J] 



th 



\PjM^- ■Pj,E^- ■Pj,S^- ■Pj,W^-' 



Suppose that position i belongs to player West and in the i position there is a black card. In this case: 

k ^ 

\D' : {b,f,ps,p4,--- ,p-k)=uj\^ bw -{s^ 1)1 -s^^Yl 

fJlPj,Nl-Pj,El-Pj,s^--Pj,wl 

With these connections and with (|3| we get the following: 

\D' : (5, f,p3,p4, • • • ,pk) = uj,m the i"^ position there is a black card| 

= \D' : (6, f,p3,p4, ■ ■ ■ ,pk) = u]\ ■ f — Ijposition i belongs to player North in the dealing method} 

v s 

bs 

H liposition i belongs to player East in the dealing method} 

s 

bs . . 

H Imposition i belongs to player South in the dealing method} 

s 

+ ^^^1 {position i belongs to player West in the dealing method} 
s 

So the most important question is which position belongs to player North, player East, player South and 
player West, and we have to analyze these positions by groups. 

First we analyze those positions which belong to player North in the dealing method. There are s such 
positions. Let position i be one of these s positions. 

E([/(i)|in the i*'' position there is a black card in permutation D' , (b,f,p3,p4, ■ ■ ■ ,pk) = 

= P (first position is non-red and non-black | in the i"^ position there is a black card in permutation D' , 

{b,r,P3,P4,- ■ ■ ,Pk) 

-|-P(second position is non-red and non-black|in the i*'' position there is a black card in permutation D' , 

{b,r,P3,P4,- ■ ■ ,Pk) ^(^) 

+ • • -^-P (position (i — 1) is non-red and non-black | in the i'^ position there is a black card in permutation D' , 

{b,r,P3,P4,- ■ ■ ,Pk) = 

— P(position (i -t- 1) is non-red and non-black|in the i^^ position there is a black card in permutation D' , 

{b,r,P3,P4,- ■ ■ ,Pk) ^(^) 

— P(position (i -I- 2) is non-red and non-black|in the i^^ position there is a black card in permutation D' , 

{b,r,P3,p4,- ■ ■ ,Pk) = 

— • • — P((position (4s) is non-red and non-black|in the i^^ position there is a black card in permutation D' , 
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(5,r,p3,P4, • • • ,Pk) = w) - 2i. 
We can compute one of these probabilities in the following way: 

P (position i is non-red and non-black | in the i'^ position there is a black card in permutation D' , 

{b,r,P3,P4,- ■ ■ ,Pk) ^(^) 

1 {position j belongs to player North in the dealing method} 



s - 1 
, Pe 



1 {position j belongs to player East in the dealing method} 
Ps 

H l{position j belongs to player South in the dealing method} 

s 

-I- ^'^l {position j belongs to player West in the dealing method}. 
s 

When calculating the conditional expectations this probability has a positive sign if j < i, and a negative 
sign if i < j. This is true for every position j ^ i, and i that belongs to player North. 

If we analyze just those positions j which belong to player East, this probability will have a positive sign 
whenever j < i, and a negative sign whenever i < j. This is what we called Z{E, N). 

If we analyze just those positions j which belong to player South, this probability will have a positive 
sign whenever j < i, and a negative sign whenever i < j. This is what we called Z{S, N). 

If we analyze just those positions j which belong to player West, this probability will have a positive sign 
whenever j < i, and a negative sign whenever i < j. This is what we called Z{W,N). 

The value Z{N, N) is trivially for each dealing method. 

Using these results we get: 



(X! 1-^' • (^''''-P3,P4, • ■ ■ ,Pk) = w, 



i=l 



in the i*'' position there is a black card in permutation D' and position i belongs to player North | 
E(J7(z)|in the z*'' position there is a black card in permutation D' and position i belongs to player North, 

ib,r,p3,pi, ■ ■ ■ ,pk) = 

inw?; I f bN ( PE ■ Z{E, N)+ps- Z{S, N)+pw Z(W, N) ' 

= \D : {b,r,p3,p4, ■ ■ ■ ,pk) = — ( 2 A^p 



So here we analized just those positions which belong to player North. 
If we analyze every position, we get the following: (see also (|3|) 

^ Z{D',B,R) 

= sliU — i-^ -) ((4.+l)5+ ^ (^^ ■ ^) + ■ ^) + ■ ^(^^ ^) „2 V N, 

^fL\Pt,Nl ■ Pi,El ■ Pi,sl ■ P^,wy ^ S V S " 

^ bE ( PN ■ Z{N,E) +PS ■ Z{S,E) +PW ■ Z{W,E) ^ 

* p=i 
^ bs ^ PN ■ Z{N,S) +PE - Z{E,S) +PW ■ Z{W,S) 

, bw ^N-Z{N,W)+pE-Z{E,W)+ps-Z{S,W) ^ ^ 



s 

so the proof is complete. 



□ 
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4 The case of more players 



Consider the following case: there are I players and, in the deck, there are £ • s cards and there are k different 
card values. 



The analogue of Q now reads as: 



J2d':{pi,P2,--- ,p-k)=uj Z{D', Pi, Pi+l) 

~{ ■ Pi+i 



+ 0(a-2). 



where pi,p2, - • • ,Pk are the number of cards coloured Pi, P2, ■ ■ ■ , Pk in the deck. 

The main question is how we can calculate Ed' (pi p2 ••• Pk)=io Z{D' , X,Y), where X and Y are two 
different types. Let Xi be the number of X coloured cards which player i receives after the dealing method 
and X be the number of cards which have value X. Let p^ j be the number of coloured cards which player 
j receives after the dealing method. 



Theorem 4.1. 



E ZiD',X,Y) 

D':(pi,p2,--- ,Pk)=uj 

1 / ((Eo=iPoj) - - Vj) ■ z{i, i) 



= \D \ [pi,p2,- ■■ ,Pk) = t^\[{ts + l)x+ -2^Xi\^ 2^ 

1=1 t=i 

= \D':ip-i,p-2,---,Pk)=u.\{±Xi{is + l-2^+ ± ((^^^Po^^)-^-y^y^(^^^ 

i=l j = l,j:^i 

where ii is the first position which goes to player i in the dealing method, 12 is the second position which goes 
to player i in the dealing method, etc. 

Proof. One can prove it in the same way as in the previous case, when there were only 4 players. □ 

Theorem 4.2. The coefficient of a^^ is exactly s times larger in the ordered dealing than in the cyclic 
dealing for every possible k and £ values, and if s is an even number then the coefficient is in the back 
and forth dealing, if s is an odd number then the coefficient is exactly s times smaller in the back and forth 
dealing than in the cyclic dealing for every possible k and I values. 

Proof. In the ordered dealing, without loss of generality, suppose that player j receives all his cards before 
player i receives his first card. In this case in the representing sequence there are s'^ j — i pairs and i — j 
pairs. So Z{j, i) = in this case. 

In the representing sequence of the cyclic dealing the first position which belongs to player j stands in 
j — i pair with s positions which belong to player i, the second position which belongs to player j stands in 
j — i pair with s — 1 positions which belong to player i, the last position which belongs to player j stands in 
j — i pair with 1 position which belongs to player i; if j receives his first card before player i. Similarly, the 
first position which belongs to player i stands in i — j pair with s — 1 positions which belong to player j, the 
second position which belongs to player i stands in i — j pair with s — 2 positions which belong to player j, 
the last position which belongs to player i stands in i — j pair with positions which belong to player j; if 
j receives his first card before player i. So Z{j,i) = s in this case. 

If s is an even number then the representing sequence of the back and forth dealing is symmetric, so 
Z{j,i)^0. 

If s is an odd number then in the representing sequence of the back and forth dealing let the first £s — £ 
positions be the first group, the last £ positions be the second group. The first group is symmetric, so within 
the first group Z{j, i) — 0. In the first group the s — 1 positions which belong to player j stand in j — i pair 
with the position which belongs to player i in the second group. In the first group the s — 1 positions which 
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belong to player / stand in i — j pair with the position which belongs to player j in the second group. If 
player j receives his first card before player i then in the second group the position which belongs to player 
j stands in j — i pair with the position which belongs to player i in the second group. So Z{j, i) = 1 in this 
case. 

So we proved the following: 
Ordered dealing: Z{j,i) = ■ Ij^i. 

Cyclic dealing: Z{j,i) = s ■ Ij^i, 

Back and forth: Z{j, i) = Ij^i, if s is an odd number. 
Back and forth: Z{j, i) = 0, if s is an even number. 
Here: 

1, player j receives his first card before player i 
— 1, player i receives his first card before player j. 

This is true for each pair of players, therefore the same holds for X]j=i j^i ^^^°=^^''-'^~2'~^'^'^^'''^^ : 
this is s times bigger for the ordered dealing than for the cyclic deahng, and this sum is s times bigger for 
the cyclic dealing than for the back and forth dealing with odd s values, while the sum is zero for the back 
and forth dealing with even s values. 

Now we analyze 2 ^*=^ ^' . 

Ordered dealing: 

In the representing sequence the position (?' — l),s + 1, (i — l)s + 2, • • • , is belong to the i"' player, so we 
have to sum up these positions when we compute J2t=i ^t- 

^^{{i-l)s + l+is)s 
s 2s 

Cyclic dealing: 

In the representing sequence the position i, 1 + i,2i + i, (s — + i belong to the i*'' player, so we have 
to sum up these positions when we compute J2t=i ^t- 

2^=2Mi^=2i + .(.-l). 
s 2s 

Back and forth dealing: 

In the representing sequence the position i, 2£ — (i — 1), 2£ + i, 4^ — (i — 1), 4£ + i, 6^ — (i — 1), • • • , £s — (i — 1) 
belong to the i*'' player if s is an even number, so we have to sum up these positions when we compute 

if s is an even number. 

In the representing sequence the position i, 2£ — (i — 1), 2£ + i, 4^ — — 1), 4£ + i, 6£ — (i — 1), • • • ,£s—{i—i) 
belong to the i*^ player if s is an odd number, so we have to sum up these positions when we compute Ylt=i 

,ELili = V(i + ^(. i))(.-i)^ ^ 1 

s s V 2 J s 



if s is an odd number. 

Now we can easi' 
Ordered dealing: 



Now we can easily analyze Is + 1 — 2 ^*-^ ^' . 



1 + 1-2 ^*=^^* =s{£-2i + l). 



Cylic dealing: 

£s + l- 2^*=!^ = _ 2i + 1). 



s 



11 



Back and forth dealing: 

s 

if s is an even number. ^ 

es + 1- 2^*^ ^ _ 2i + 1), 
s s 

if s is an odd number. 

Thus we see that these terms also differ by factors of s when comparing the ordered, cyclic, and back 
and forth dealing methods (odd s values), while this term is also for the back and forth dealing if s is an 
even number. 

We proved the claim for each term, so it is true for J2d'-(pi p2 ■■■ Pk)=^^ X, Y) and therefore it is also 

true for the coefficient of a~^, so we proved the statement. 

□ 



5 The case of arbitrary initial deck 

Now, we suppose that the initial deck is arbitrary. In this case the variation distance is the following: 



\3 = l 



E 

a<b 



WiD, a, b) ED':fe,p„... ,pO=^ Z{D', a, h) 



ria ■ rib 



+ 0{a-^). 



In this formula the only term which depends on the dealing method is pa ••• Pk)=i^ -^(^'i the 

previous section we proved for these terms the following: these terms differ by factors of s when comparing 
the ordered, cyclic, and back and forth dealing with odd s values, while this term is also for the back and 
forth dealing if s is an even number. We can easily see that these terms do not depend on the initial deck, 
therefore for every initial deck the coefficient of a^^ is s times bigger for the ordered dealing than for the 
cyclic dealing and s times bigger for the cyclic dealing than for the back and forth dealing with odd s values, 
while the coefficient of a"^ is zero for the back and forth dealing with even s values. 

6 The case of two or three different types of cards 

The purpose of this section is to gain some quantitative insight on how the leading term of the variation 
distance behaves in the case of repetitive cards. 



6.1 Two different types of cards 

In this section we analyze those cases when there are only two types of cards: black and red, 52 cards in the 
deck and four players. Using our results we know that the coefficient of is exactly 13 times smaller in 
the cyclic dealing method than in the ordered dealing method, and 13 times smaller in the back and forth 
dealing method than in the cyclic dealing method. Let h be the number of black cards in the deck. Using 
the computer and our formulas we are able to compute the coefficient of a^^ for any possible value b. 

the coefficient of in mI) will be the following: 



With using Theorem 3.1 



13 



6(52 



E 



26jv 
13 




13 , 13 

p—i p—i 

In this formula there are four sums: the sum of positions i where f{i) = N, the sum of positions i where 
f{i) = E, the sum of positions i where f{i) = S and the sum of positions i where f{i) — W. In the ordered 
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dealing the values of these sums are the following: 91; 260; 429; 598. In the cyclic dealing the values of these 
sums are the following: 325; 338; 351; 364. In the back and forth dealing the values of these sums are the 
following: 343; 344; 345; 346. We see that the distances of the different values of the sums are the smallest 
in the back and forth dealing and the largest in the ordered dealing in these methods. In this sense the best 
dealing of these methods is the back and forth dealing, and the worst method is the ordered dealing. We can 
check that in an unfeasible dealing method in which the sums are 344,5; 344,5; 344,5; 344,5 the coefficient of 
would be 0. We think that dealing method A is better than dealing method B if the following value is 
smaller in A than in B: \ J^ll, Np - 344, 5| + | ^pli Ep - 344, 5| + | ^pli Sp - 344, 5| + | i Wp - 344, 5|. 
Our conjecture is that the best dealing method is a dealing method in which two sums are 344 and two sums 
are 345. As an example, let g be the following: 

5(2) = .9(8) = .g(9) = .g(16) = .g(17) = g{24) = g{25) - g{32) = ^(33) = g(40) 
- .g(41) - .g(48) - .g(49) = TV, 

5(3) = 5(7) = 5(10) = 5(15) = 5(18) = 5(23) = 5(26) = 5(31) = 5(34) = 5(39) 
= 5(42) = 5(47) = 5(50) = E, 

5(1) = ff(6) = 5(12) = 5(14) = 5(19) = 5(22) = 5(27) = 5(30) - 5(35) = 5(38) 
= 5(43) = 5(46) = 5(51) = S, 

5(4) = 5(5) - 5(11) = .9(13) - 5(20) - 5(21) - 5(28) - 5(29) - 5(36) - 5(37) 
= g(44) = 5(45) = 5(52) = W. 
So the representing sequence is: 

SNEWWSENNEWSWSENNESWWSENNESWWSEN 

NESWWSENNESWWSENNESW. 

If we sketch a graph we can see that this dealing method has a smaller coefficient than the back and forth 
dealing for each value b. 

If we compare these coefficients and the coefficients of the back and forth dealing, we can see that these 
coefficients are circa twice smaller than the coefficients of the back and forth dealing. So for large a we can 
save circa 1 riffle shuffle if we use this dealing method instead of the back and forth dealing. We note that 
there are other dealing methods which have the same coefficient of a^^ for each value b, but our conjecture 
is that we cannot construct a better dealing method. 

6.2 Three different types of cards 

Now we suppose that there are three different types of cards in the deck: black, red and green. The number 
of cards is 52 and there are four players. By our results we know that the coefficient of is exactly 13 
times smaller in the cyclic dealing method than in the ordered dealing method, and 13 times smaller in the 
back and forth dealing method than in the cyclic dealing method. Let b be the number of black cards, r be 
the number of red cards, 5 be the number of green cards in the deck. If we know b and r we also know 5, 
because 5 = 52 — 6 — r. Using the computer and our formula we are able to compute the coefficient of 
for any possible b, r values. Figure [2] shows the coefficient of a^^ for each possible b, r values in the back and 
forth dealing method. 

A main question is whether the former dealing method 5 in which two sums are 344 and two sums are 
345 is better than the back and forth dealing in this case. The answer is: not for every configurations. For 
configuration: & = l,r = 1 the coefficient of a^^ is for the black and forth dealing and for dealing 
method 5. So we can see that for this configuration the back and forth dealing is better than dealing method 
5. Figure [3] shows the coefficient of a^^ for each possible b,r values in dealing method 5. In most cases 
dealing method 5 has a smaller coefficient than the back and forth dealing, but there are some configurations 
when the back and forth dealing is better than dealing method 5. This proves that the conjecture what we 
drew up for two different types is false for this case. Hence this conjecture is not true for arbitrary different 
types of cards. 
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Figure 1: The coefficient of for each possible value b in the back and forth and in our conjectured best 
dealing method. The horizontal axis marks the value b and the vertical axis marks the coefficient of a~^. 
Square plots our conjectured best dealing method and circle plots the back and forth dealing method. For 
better illustration we did not represent the points b — 1 and 51 for the ordered dealing. 
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Figure 2: The coefficient of a ^ for each possible value b, r in the back and forth dealing method. 




